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ABSTRACT 

The chiral glass behaviour of the nearest-neighbour random-bond XY spin 
glass in four dimensions is studied by Monte Carlo simulations. A chiral 
glass transition at T cg = 0.90 ±0.05 is found by a finite-size scaling analysis 
of the results. The associated chiral correlation-length exponent is esti- 
mated to be v cg = 0.6 ± 0.1 and rj cg ~ 0.25. The values for the chiral 
critical temperature and the exponents are very similar to those reported 
recently for the spin glass transition in this model. The results strongly 
suggest a simultaneous ordering of spin and chirality in four dimensions. 
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1. Introduction 

Although it has been known since Villain [1] that frustrated vector spin 
systems possess both reflectional and rotational symmetries, it's only fairly 
recently that the chiral glass behaviour of vector spin glasses has been stud- 
ied [2-9]. Whereas the continuous rotational symmetry is associated with 
the spins, it is the discrete twofold Ising-like reflectional symmetry which is 
associated with chirality. Much of the recent work has been motivated by 
the suggestion that the chiral glass transition in vector spin glasses belongs 
to the same universality class as the Ising spin glass transition. This would 
imply that the spin glass transition observed in experiments may, in fact, 
be 'chirality driven' [7]. Consequently, one would have a chiral glass phase 
with broken reflectional symmetry but preserved rotational symmetry. Re- 
cent numerical work by Kawamura [4,7,8] in three dimensions would appear 
to be consistent with this picture. 

There is now convincing evidence that both XY [3-6,10] and Heisenberg 
[6-8,11] spin glasses exhibit conventional spin glass ordering only at zero 
temperature for d = 2 and 3. Domain- wall renormalization-group studies 
suggest chiral ordering for both XY [5] and Heisenberg [7] spin glasses in 
two dimensions also at zero temperature only. However, the values of the 
chiral and spin glass correlation-length exponents have been found to differ 
in both cases. Further evidence for a decoupling of the chiral and phase 
variables on long length scales has come from various Monte Carlo studies 
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in 2d [2,3], including very recent work using a vortex representation [9]. 
In 3d the results point to a finite-temperature chiral glass transition both 
for XY [4-6] and Heisenberg [6-8] spin glasses (although the case is far from 
convincing for the latter) . So it would appear that chiralities and spins have 
markedly different behaviour in both d = 2 and 3. 

Very recent Monte Carlo simulations in four dimensions suggest a finite- 
temperature spin glass transition for vector spin glasses [12,13]. This would 
imply that the lower critical dimensionality for vector spin glasses is less 
than four. It is clearly of interest to see whether the difference in the 
behaviour of chiralities and spins extends into higher dimensions. 

In this paper we present the results from Monte Carlo simulations of 
the chiral behaviour of the four dimensional random-bond XY spin glass. 
Using finite-size scaling and our earlier results for the spin glass transition 
[12], we shall find evidence for a simultaneous ordering of spin and chirality 
in four dimensions. 

In section 2 we define the model and review the finite-size scaling technique 
used to analyse the data. Section 3 gives details about the simulations and 
the results are presented and discussed in section 4. The conclusion is given 
in section 5. 
2. The model 

The Hamiltonian for the model is given by 

H = ~Y1 Jijcoatfi-ej), (1) 

<i,j> 
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where < Q{ < 2n for all planar spins i and the summation runs over 
all nearest-neighbour pairs on a four dimensional hypercubic lattice with 
L x L x L x L spins (L = 2, 4 and 6) with full periodic boundary con- 
ditions. The interactions, Jij, are independent random variables selected 
from a binary ±1 distribution. Throughout, the temperature is given in 
units of the nearest-neighbour interaction. We recently reported the results 
concerning the spin glass transition for this model [12]. In this publication 
we concentrate on the chiral glass transition. 

The local chirality, K a , at a plaquette a consisting of four spins is 
defined by the scalar [2] 



where the summation is performed over a directed (clockwise) closed con- 
tour along the four sides of the plaquette. Clearly, whereas n a = for an 
isolated unfrustrated plaquette, it's restricted to the values ±1 for frustrated 
ones. Hence, chirality act as a 'continuous' Ising-like quantity. However, as 
it's magnitude fluctuates with the temperature to a certain extent, we work 
with a root-mean-square value given by [2] 



where < ... >t denotes a thermal average, [...]j indicates an average over 
the disorder and Nd, the total number of plaquettes for our four dimensional 




(2) 



a 




(3) 



lattice , is equal to 6L 4 . 
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At a chiral glass transition one expects the chiral glass susceptibility, Xcg, 
to diverge, where 



Xcg 



j 



a,/3 



eg •> 



(4) 



and here the summation is with respect to all plaquettes a and (3. 

(2) 

The chiral glass order parameter, q cg \ defined above can be written in terms 
of the overlap between two replicas 1 and 2, viz 



where 



Qcf = [< Q 2 >t]j, 



1 1 2 



(5) 



(6) 



a 



In order to improve the analogy with Ising spins, we follow Kawamura [2] 
and work in the simulations with the reduced chiral glass susceptibility 



Xcg — Xcg I K ■ 



(7) 



Another key quantity studied in the simulations is the dimensionless Binder 
parameter defined by [14] 



and here 



9c g (L,T) = ± 



,( 4 ) - \s „ 4 



(4) 

Teg 



(8) 



teg 



[< q >t]j- 



(9) 
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The Binder parameter is expected to scale as [14] 

g cg (L,T)=g cg (L 1 ^(T-T cg )), (10) 

where T cg and v cg are the chiral glass critical temperature and correlation- 
length exponent, respectively. The value of T cg can be located by using the 
fact that the scaling function, g cg , and, hence, also g cg are independent of L 
at the transition temperature. The chiral glass correlation-length exponent 
then follows from a one-parameter scaling fit of the data for g cg . 
The finite-size scaling form for the reduced chiral glass susceptibility is given 

by 

Xc 9 (L,T) =L 2 -^ cg (L 1 ^(T-T cg )). (11) 

Here rj cg is the chiral critical-point decay exponent and x cg the scaling 
function. 

The finite-size scaling form for the analogous spin glass quantities can 
be found in [12]. 
3. Simulations 

We now discuss the computer simulations (Jain [12] should be consulted 
for further technical details). In order to ensure that equilibrium has been 
achieved in the simulations, we use the technique of Bhatt and Young [14] 
whereby chiral glass correlations (qig and qig ) computed from two replicas 
at the same time are required to agree with those from one replica at two 
different times. 
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At low temperatures the chiral degrees of freedom were more difficult to 
equilibrate than their spin counterparts, requiring approximately twice as 
many Monte Carlo steps. As a result, the lowest temperatures studied were 
T = 0.8(L = 6), 0.7 (L = 4) and 0.3(L = 2). We took disorder averages over 
many pairs of independent samples: 100 ~ 120(1/ = 6), 100 ~ 250(1/ = 4) 
and 250 ~ 700(1/ = 2) for each temperature. Most of the computational 
time (about 200 hours of CPU time on a Cray J932) was taken up by 
the simulations for L = 6. (The processors of the Cray J932 are roughly 
only half the speed of those of the Cray YMP.) It's estimated that an 
additional 1000 hours of CPU time would be required to obtain reliable 
data on a larger lattice such as L = 8. All computational studies performed 
to-date, including the present one, on four-dimensional spin glasses have 
been restricted to small lattices with L < 6 [12-14]. 

4. Results 

In this section we present the results. 

Figure 1 shows a plot of the reduced chiral glass susceptibility against tem- 
perature for the three different lattices considered in this work. The statis- 
tical error-bars have been estimated from the sample-to-sample fluctuations 
and, in most cases, are smaller than the size of the data points. 

We see that the value of Xcg remains quite low until T « 1.0 and then 
increases rapidly with the system size. In fact, whereas for low temperatures 
Xcg is an increasing function of L, for higher temperatures it appears to 
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actually decrease with increasing L. This suggests that the scaling regime 
is probably quite narrow. Very similar behaviour was found recently by 
Kawamura in 3d [4]. 

In figure 2(a) we show a plot of g cg against the temperature for 0.4 < T < 
1.5. The behaviour of the Binder parameter for high temperatures T > 1.0 
indicates a disordered chiral phase. For T < 1.0 there is a sharp increase 
in g cgi indicating a build up of chiral correlations. The increase is more 
noticeable for the larger lattices. We notice that for T > 1.0 g cg assumes 
neagtive values. A negative Binder parameter has also been seen for both 
XY [4] and Heisenberg [8] spin glasses in 3d. This feature would appear to 
be a consequence of the fact that n a = on unfrustrated plaquettes even 
in the ordered state. 

Figure 2(b) displays on a much expanded scale the data in figure 2(a) 
in the vicinity of the important temperature region (0.8 < T < 1.0). We see 
that the curves appear to intersect at T cg « 0.9. Furthermore, below T cg 
the values corresponding to L = 6 are consistently above those for L = 4 
and the curves clearly splay out. We estimate the chiral glass transition 
temperature to be T cg = 0.90 ± 0.05. This value is very close to the spin 
glass transition temperature (T sg = 0.95 ± 0.15) we reported recently for 
this model [12]. Our conclusion is based on the results for small lattices and 
it is highly desirable to obtain additional data for larger lattices to confirm 
our findings. 
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Further evidence for a finite-temperature transition comes from a scaling 
plot of the Binder parameter. By varying the value of v cg and considering 
values 0.85 < T cg < 0.95, we estimate the chiral correlation-length exponent 
to be v cg = 0.6 ± 0.1. It should be noted that the error-bar quoted here is 
simply an estimate that demarcates the range of values beyond which the 
data do not scale well. Figure 3 shows the data for g cg against L}^ Vc a{T — 
T cg ) with T cg = 0.90 and v cg = 0.6. This scaling plot is far better than the 
corresponding plot for g sg (see figure 4 in [12]) and the data (including those 
for L = 2) would appear to scale particularly well near T cg . Once again, 
the value of v cg would appear to be very similar to that of v sg = 0.70 ±0.10 
[12]. 

We have appreciable uncertainty in the values of both T cg and v cg . Fur- 
thermore, the increase in the reduced chiral glass susceptibility for T < 1.0 
is extremely sharp. As a consequence, the critical region is very nar- 
row. To obtain an estimate for the decay exponent, j] cgi we tried vari- 
ous different possible values of T cg and v cg . The best such scaling plot is 
shown in figure 4 where we display Xcg/L 2-7109 against L l ' Vc s (T — T cg ) with 
T cg = 0.85, v cg = 0.8 and r/ cg = 0.25. 

Clearly, the Xcg data for L = 2 do not scale all that well, especially for 
the higher temperatures. Nevertheless, we note that our estimate for rj cg is 
not incompatible with the value for rj sg found earlier [12]. 
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5. Conclusion 

To conclude, we have presented the results of a Monte Carlo simulation of 
the chiral glass behaviour of the four dimensional random-bond XY spin 
glass on small lattices (L < 6). By means of a finite-size scaling analysis of 
the data, we have estimated both the chiral transition temperature and the 
critical exponents. The chiral glass values are very similar to their spin glass 
counterparts for this model (T cg « T sg ,v cg « v sg and rj cg « rj sg ). Hence, 
our results strongly suggest a simultaneous ordering of spin and chirality in 
four dimensions. 

Further data on larger lattices (L > 6) is required to confirm both the 
transition temperature and the critical exponents. 
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FIGURE CAPTIONS 



Figure 1 

A plot of the reduced chiral glass susceptibility, Xcg, against the temperature 
for L = 2, 4 and 6 (see equations (4) and (7) in the text). The lines are just 
to guide the eye. 

Figure 2(a) 

A plot of the Binder parameter g cg against the temperature for L = 2, 4 
and 6. 

Figure 2(b) 

Same as figure 2(a) but on a much expanded scale around the interesting 
temperature region (T « 0.9). The lines are just to guide the eye. 

Figure 3 

A scaling plot of the Binder parameter g cg versus I/ 1 /°- 6 (T — 0.9). The line 
is just a guide to the eye. 
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Figure 4 

A scaling plot of Xcg/L 2 ^ against L 1 ^^ (T - T cg ) with T cg = 0.85, v cg = 
0.8 and rj cg = 0.25 (see equation (11) and the text). The line is just a guide 
to the eye. 
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